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Abstract 

Microgrids have been identified as key components of modern electrical systems to facilitate the integration of renewable 
distributed generation units. Their analysis and controller design requires the development of advanced (typically 
model-based) techniques naturally posing an interesting challenge to the control community. Although there are widely 
accepted reduced order models to describe the dynamic behavior of microgrids, they are typically presented without 
details about the reduction procedure—hampering the understanding of the physical phenomena behind them. Preceded 
by an introduction to basic notions and definitions in power systems, the present survey reviews key characteristics and 
main components of a microgrid. We introduce the reader to the basic functionality of DC/AC inverters, as well as 
to standard operating modes and control schemes of inverter-interfaced power sources in microgrid applications. Based 
on this exposition and starting from fundamental physics, we present detailed dynamical models of the main microgrid 
components. Furthermore, we clearly state the underlying assumptions which lead to the standard reduced model with 
inverters represented by controllable voltage sources, as well as static network and load representations, hence, providing 
a complete modular model derivation of a three-phase inverter-based microgrid. 

Keywords: Microgrid modeling, microgrid analysis, smart grid applications, inverters 


1. Introduction 

1.1. Motivation 

It is a widely accepted fact that fossil-fueled thermal 
power generation highly contributes to greenhouse gas emis¬ 
sions [Ml [SUEZ]- In addition, a growing stream of scien¬ 
tific results eh mo mi has substantiated claims that these 
emissions are a key driver for climate change and global 
warming. As a consequence, many countries have agreed 
to reduce their greenhouse gas emissions. 

Apart from a reduction of energy consumption, e.g ., 
through an increase in efficiency, one possibility to reduce 
greenhouse gas emissions is to shift the energy production 
from fossil-fueled plants towards renewable sources [66] [67] 
m ■ Therefore, the worldwide use of renewable energies 
has increased significantly in recent years m- 
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Unlike fossil-fueled thermal power plants, the major¬ 
ity of renewable power plants are relatively small in terms 
of their generation power. An important consequence of 
this smaller size is that most of them are connected to the 
low voltage (LV) and medium voltage (MV) levels. Such 
generation units are commonly denoted as distributed gen¬ 
eration (DG) units pQ. In addition, most renewable DG 
units are interfaced to the network via DC/AC inverters. 
The physical characteristics of such power electronic de¬ 
vices largely differ from the characteristics of synchronous 
generators (SGs), which are the standard generating units 
in existing power systems. Hence, different control and 
operation strategies are needed in networks with a large 
amount of renewable DG units [40. 104[ llOll . 

1.2. The microgrid concept 

One potential solution to facilitate the integration of 
large shares of renewable DG units are microgrids EH HE, 
0DJ EH HI [M] ■ A microgrid gathers a combination of gen¬ 
eration units, loads and energy storage elements at distri¬ 
bution or sub-transmission level into a locally controllable 
system, which can be operated either in grid-connected 
mode or in islanded mode, i.e., in a completely isolated 
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manner from the main transmission system. The micro¬ 
grid concept has been identified as a key component in 
future electrical networks [H5J [321@21 @21 • 

Many new control problems arise for this type of net¬ 
works. Their satisfactory solution requires the develop¬ 
ment of advanced model-based controller design techniques 
that often go beyond the classical linearization-based nested- 
loop proportional-integral (PI) schemes. This situation 
has, naturally, attracted the attention of the control com¬ 
munity as it is confronted with some new challenging con¬ 
trol problems of great practical interest. 

It is clear that to carry out this task it is necessary 
to develop a procedure for assembling mathematical mod¬ 
els of a microgrid that reliably capture the fundamental 
aspects of the problem. Such models have been devel¬ 
oped by the power systems and electronics communities 
and their pertinence has been widely validated in simu¬ 
lations and applications |2B E31 ED EU- However, these 
are reduced or simplified, i.e.. linearized, models that are 
typically presented without any reference to the reduction 
procedure—hampering the understanding of the physical 
phenomena behind them. 

1.3. Existing literature 

For the purposes of this survey, previous work on mi¬ 
crogrid modeling can be broadly categorized into two classes. 
The first class focusses on modeling and control of inverter- 
interfaced DG units in microgrid applications, but the 
model derivation is restricted to individual DG units and 
the current and power flows between different units are 

not considered explicitly [521IHT11551 [TT21IM1I5T1ITU1 ITT] ■ 
The second class discusses models of microgrids includ¬ 
ing electrical network interactions, but the model deriva¬ 
tion is based on linearization (z.e., the so-called small- 
signal model) [5411811174] , Furthermore, this class of mod¬ 
eling is often tied to specific network control schemes, 
such as droop control [20] [HI] or to specific test networks 
[S3 EH [7H [73]. Building on this previous work and in a 
survey-like manner, the present paper brings both afore¬ 
mentioned classes together to formulate a generic modular 
model of a microgrid. 

Going beyond a mere review of existing microgrid mod¬ 
els, we employ model reduction via a time-scale separation 
together with the subsequent derivation of the well-known 
power flow equations, which is a standard procedure in 
SG-based networks jl05] . A similar approach has also 
been employed for microgrids in [H [72 111 IM]- How¬ 
ever, neither reference provides a detailed model derivation 
for inverter-interfaced units. Also, the analysis in ciin] 
is restricted to an AC microgrid consisting of two invert¬ 
ers connected via a resistive-inductive line and two local 
resistive-inductive loads, while the modeling procedure of 
the present paper applies to networks with generic topo¬ 
logy and arbitrary number of units. 


1.4. About the survey 

The present survey is an attempt to provide a guide¬ 
line for control engineers attracted by this fundamental 
application for Smart Grids to assess the importance of 
the main dynamical components of a three-phase inverter- 
based microgrid as well as the validity of different models 
used in the power literature. To this end, we at first re¬ 
view some fundamental concepts and definitions in power 
systems, including a survey on the notion of instanta¬ 
neous power. Subsequently, we introduce the reader to 
the microgrid concept and discuss its main components. 
We illustrate that inverter-interfaced units are the main 
new elements in future power networks, detail the ba¬ 
sic functionality of inverters and review the most com¬ 
mon operation modes of inverter-interfaced units together 
with their corresponding control schemes. This paves the 
path for—starting from fundamental physics—presenting 
detailed dynamical models of the individual microgrid com¬ 
ponents. Subsequently, we clearly state the underlying 
assumptions which lead to the standard reduced model 
with inverters represented by controllable voltage sources, 
as well as static network and load representations. This 
reduced model is used in most of the available work on 
microgrid control design and analysis [96l[90U2l[76U29] . 

We focus on purely inverter-based networks, since in¬ 
verter-interfaced units are the main new elements in mi¬ 
crogrids compared to traditional power systems. However, 
we remark that the employed modeling and model reduc¬ 
tion techniques can equivalently be applied to standard 
bulk power system models as well as to power systems 
with mixed generation pool. For modeling of traditional 
electro-mechanical SG-based units, the reader is referred 
to standard textbooks on power systems [521 [521 [5], 

The main contributions of the present survey paper are 
summarized as follows. 

• Provide a detailed comprehensive model derivation 
of a microgrid based on fundamental physics and 
combined with detailed reviews of the microgrid con¬ 
cept, its components and their main operation modes. 

• Answer the question, when an inverter can be mod¬ 
eled as a controllable AC voltage source and depict 
the necessary underlying model assumptions. 

• Show that the usual power flow equations can be 
obtained from a network with dynamic line models 
via a suitable coordinate transformation (called dq- 
transformation) together with a singular perturba¬ 
tion argument. 

• By combining the two latter contributions, recover 
the reduced-order microgrid model currently widely 
used in the literature. 

We emphasize that the aim of the present survey is not 
to give an overarching justification for the final (simpli¬ 
fied) model, but to provide a comprehensive overview of 
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the modeling procedure for main microgrid components 
together with their dynamics, as well as of the main nec¬ 
essary assumptions, which allow the reduction of model 
complexity. Which of the presented models (if any) is ap¬ 
propriate for a specific control design and analysis cannot 
be established in general, but has to be decided by the 
user. Any model used in simulation and analysis necessar¬ 
ily involves certain assumptions. Therefore, it is of great 
importance that the user is aware of the pertinence of the 
employed model to appropriately assess the implications 
of a model-based analysis. 

The remainder of this survey is structured as follows. 
Basic preliminaries, such as definitions and transforma¬ 
tions in power systems are given in Section [2] The micro¬ 
grid concept is reviewed in Section[3] A detailed dynamical 
model of a microgrid is derived in Section]?] In particular, 
common operation modes of inverter-interfaced units are 
discussed therein. The model reduction yielding models of 
inverters as AC voltage sources and phasorial power flow 
equations is conducted in Section [5] The survey paper 
is wrapped-up with some conclusions and a discussion of 
topics of future research in Section [hj 

Notation. We define the sets R>o := {x £ R|x > 0}, 
R>o := {x £ R|x > 0} and S := [0, 2ir). For a set V, 
let |V| denote its cardinality and [V] fc denote the set of all 
subsets of V that contain k elements. For a set of, possibly 
unordered, positive natural numbers V = {l, k ,..., n}, the 
short-hand i ~ V denotes i = l,k,... ,n. Given a positive 
integer n, we use 0„ to denote the vector of all zeros, 1„ 
the vector with all ones and I„ the n x n identity matrix. 
Let x = col(xi,..., x n ) £ C n denote a column vector with 
entries Xi £ C. Whenever clear from the context, we simply 
write x = col(xj) £ C n . The 2-norm of a vector x £ C n 
is denoted by ||x ||2 = -y/|xi| 1 2 + ... + \x n \ 2 . Let diag(ai) £ 
C" xn denote a diagonal matrix with entries a,; £ C. Let 
j denote the imaginary unit. The conjugate transpose of 
a matrix A £ C nx " is denoted by A*. For a function / : 
R” —► R, V/ denotes the transpose of its gradient. The 
operator ® denotes the Kronecker product. 

2. Preliminaries and basic definitions 

2.1. Symmetric AC three-phase signals 

Definition 2.1. \27^ A signal x : R>o —> R is said to be 

an AC signal if it satisfies the following conditions 

1. it is periodic with period T £ R>o, he., 

x(t) = x(t + nT), Vn £ N, Vi > 0, 

2. its arithmetic mean is zero, i.e., 

rt-\-T 

J x{r)dr = 0 Vt > 0. 


Definition 2.2. A signal x : R>o —> R 3 is said to be a 
three-phase AC signal if it is of the form 

xabc = col(x A ,x B ,x c ), 

where xa ■ R>o —t R, x B '■ R>o —► R and xc '■ R>o —► R 
are AC signals. 

A special kind of three-phase AC signals are symmetric 
AC three-phase signals, defined below. 

Definition 2.3. |?} Chapter 2] A three-phase AC signal 
x a bc '■ R>o —t R 3 is said to be symmetric if it can be de¬ 
scribed by 



X a {t) 


*-0 

if 

'm 


x b {t) 

= A(t ) 

sin (S(t) — 


_x c (t) 


sin (S(t) + ^f) 


where A : R>o —> R>o is called the amplitude and S : 
R>o —> S is called the phase angle of the signal. 

Clearly, from the preceding definition, a symmetric 
three-phase AC signal x a bc can be described completely 
by two signals: its angle S and its amplitude]]] A. 

Definition 2.4. Q Chapter 2] A three-phase AC signal 
is said to be asymmetric if it is not symmetric. 

Definition 2.5. J^<g[ Chapter 3] A three-phase AC elec¬ 
trical system is said to be symmetrically configured if a 
symmetrical feeding voltage yields a symmetrical current 
and vice versa. 

Definition 2.6. f48\ Chapter 3] A three-phase AC power 
system is said to be operated under symmetric conditions 
if it is symmetrically configured and symmetrically fed. 

Examples of symmetric and asymmetric three-phase 
AC signal]]] are given in Fig. [l] Note that the signal in 
Fig. [lb] satisfies Definition |2.3| and only differs from the 
signal in Fig.[la]iii that it possesses a time-varying periodic 
amplitude. 

Remark 2.7. The terms “balanced” and “unbalanced” 
are frequently used as synonyms of “symmetric”, respec¬ 
tively “asymmetric” in the literature HE7]. 

Remark 2.8. Three-phase electrical power systems con¬ 
sist of three main conductors in parallel. Each of these 
conductors carries an AC current. A three-phase system 
can be arranged in A- or Y-configuration, see Fig. [2j The 


1 To simplify notation the time argument of all signals is omitted, 
whenever clear from the context. The same applies to the definition 
of signals, i.e., a signal x : R>o —>• K, is defined equally as x E R, 
whenever clear from the context. 

2 In this paper only AC systems and signals are considered. There¬ 

fore, the qualifier AC is dropped from now on. 
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(a) Symmetric three-phase AC (b) Symmetric three-phase AC 
signal with constant amplitude signal with time-varying ampli¬ 
tude 


S —> R 3x3 , 
Tdqo{g) ■= 



"cos (g) 
sin(e) 

V2 

L 2 


cos(p — 

sin (g(t) - 

V2 

2 



Then, fd q o : l 3 x § -) R 3 , 


cos (g + 

sin(g + 

VS 

2 



( 2 . 1 ) 


fd q o(x{t), g(t)) = T dq0 (g(t))x(t) (2.2) 



s 

(c) Asymmetric three-phase AC 
signal with phases not shifted 
equally by If 



s 

(d) Asymmetric three-phase AC 
signal resulting of an asymmet¬ 
ric superposition of a symmetric 
signal with signals oscillating at 
higher frequencies 


Figure 1: Symmetric and asymmetric AC three-phase signals. The 
lines correspond to x a ’— x & ’- x c ’• • • ’• 



is called dqO-transformation. 


Note that the mapping (2.1) is unitary, i.e., Tj q0 = Tf qQ . 
From a geometrical point of view, the d^O-transformation 
is a concatenation of two rotational transformations, see 
m for further details. The variables in the transformed 
coordinates are often denoted by dqO- variables. 

The dqO-transformation offers various advantages when 
analyzing and working with power systems and is therefore 
widely used in applications [77} [710} 1101111121 . For exam¬ 
ple, the dgO-transformation permits, through appropriate 
choice of g, to map three-phase AC signals to constant sig¬ 
nals, i.e., to transform periodic orbits into constant equi¬ 
libria. This simplifies the control design and analysis in 
power systems, which is the main reason why the trans¬ 


formation (2.2) is introduced in the present case. In ad¬ 


dition, the transformation (2.21 exploits the fact that, in 


a power system operated under symmetric conditions, a 
three-phase signal can be represented by two quantities. 
To see this, let x a b c ■ R>o —> R 3 be a symmetric three- 
phase signal with amplitude A : R> 0 —X R>o and phase 
angle 9 : R>o —X S, as in Definition 2.3 Applying the 
mapping (2.1) with some angle g : R>o —X S to x a b c yields 


Figure 2: Standard Y- and A-configurations of three-phase AC 
power systems based on [481 Chapter 3]. Zy £ C denotes the 
impedance in Y-configuration, while Z/± E C denotes the impedance 
in A-configuration. 


XdqO 


X d 


/3 

sin(0 — g) 

X q 

X 0 

— ^dqO^Q^XaJjc — 


cos (9 — g) 
0 


latter is also called wye-configuration. Frequently, in a sys¬ 
tem with Y-configuration an additional fourth (grounded) 
neutral conductor is used to reduce transient overvoltages 
and to carry asymmetric currents m Chapter 2]. Such 
systems are typically called three-phase four-wire systems. 
Most three-phase power systems are four-wire Y-connected 
systems with grounded neutral conductor m Chapter 2]. 
However, it can be shown that, under symmetric operat¬ 
ing conditions, this fourth wire does not carry any current 
and can therefore be neglected [57} Chapter 2]. 


Hence, Xq = 0 for all t > 0. Therefore and as in this work 
only symmetric three-phase signals are considered, it is 
convenient to introduce the mapping Td q : S —> R 2x3 , 


T d q {e) ■= 


/2 rcos(p) 
V 3 [sin(p) 


cos (g 
sin(p 


M 


which, when applied to the symmetric 
Xabc defined above, yields 


cos(p + 17r) 
sin(p + 1 7 r) ’ 

(2.3) 

three-phase signal 


%dq — 


Xd 

Xq 


Tdq(o)x abc 


l~3 r S in(f? — g) 

V 2 cos(0 — g) 


2.2. The dqO-transformation 

An important coordinate transformation known as dqO- 
transformation in the literature [751173 HE [Ml 11011 HE I112| 
is introduced. 

Definition 2.9. [5l Chapter 4], Wh Chapter 11] Let x : 
R>o —X R 3 and g : R>o —> S. Consider the mapping T ^ 0 : 


In the following, Xd q are referred to as the dq-coordinates 
of x abc . Note that x abc = T dq (g) T x dq . 


Remark 2.10. There are several variants of the mapping 
(2.1) available in the literature. They may differ from the 
mapping © in the order of the rows and the sign of the 
entries in the second row of the matrix given in (2.1), see 
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urnm- However, all representations are equivalent in 
the sense that they can all be represented by T dq o as given 
in ( 2.11 by choosing an appropriate angle g and, possibly, 
rearranging the row order of the matrix T dq o- The same 
applies to the mapping Td q given in (2.3). Since—with a 
slightly different scaling factor -this transformation was 
first introduced by Robert H. Park in 1929 [75] it is also 
often called Park transformation [1311 Appendix A]. 


2.3. Instantaneous power 

Power is one of the most important quantities in con¬ 
trol, monitoring and operation of electrical networks. The 
first theoretical contributions to the definition of the power 
flows in an AC network date back to the early 20 th cen¬ 
tury. However, these first definitions are restricted to sinu¬ 
soidal steady-state conditions and based on the root-mean- 
square values of currents and voltages. As a consequence, 
these definitions of electric power are not well-suited for 
the purposes of network control under time-varying oper¬ 
ating conditions Hj. 

The extension of the definition of electrical power to 
time-varying operating conditions is called “instantaneous 
power theory” in the power system and power electronics 
community pi non. The development of this theory al¬ 
ready begun in the 1930s with the study of active and non¬ 
active components of currents and voltages |35| . Among 
others, relevant contributions are [TS1 [22] 131 11091 12311731 

Ell- 

Today, it is widely agreed by researchers and prac¬ 
titioners PM1 [731IM] that the definitions of instanta¬ 
neous power proposed in J3| and contained in [|j are well- 
suited for describing the power flows in three-phase three- 
wire systems and symmetric three-phase four-wire sys¬ 
tems. However, a proper definition of instantaneous power 
in asymmetric three-phase four-wire systems with nonzero 
neutral current and voltage is still an open (and contro¬ 
versial) field of research [2j, ffl El lOlj ■ A good overview 
of the research history on instantaneous power theory is 
given in [Toll Appendix B], 

Consider a symmetric three-phase voltage, respectively 
current, given by 



sin(0) 


sin(ip) 

Vabc = V2V 

sin(0 - f) 

? 'iabc — V%I 

sin(<^ — vy) 


_sin(6>+ 2p)_ 


sin ( 9 ? + ^f) 


(2.4) 

where 9 : M>o —> S and ip : K>o S are the phase angles 
and \f2V : K>o — > R>o, respectively \/2I : R> 0 —> M> 0 , 
the amplitudes of the respective three-phase signal. As 
applying the transformation (2.3) to 


shown in Section 2.2 


the signals given in (2.4) yields 


Vdq = 

1 _ 1 

= 

sin(0 — g) 
cos [6 — q) 

^dq = 

1 - 1 

1 _ 1 

= VZI 

sin(</5 - g) 
cos{p — g) 


(2.5) 


Based on the preceding discussion, the following defini¬ 
tions of instantaneous active, reactive and apparent power 
under symmetric, but not necessarily steady-state, condi¬ 
tions are used in this work. The definitions are based on (31 
[4], in which they are given in a/3-coordinates. For the pur¬ 
pose of the present paper, it is more convenient to equiv¬ 
alently define the instantaneous powers in dg-coordinates. 


Definition 2.11. Let Vd q {t) and i dq (t) be given by (2.5). 
The instantaneous three-phase active power is defined as 


P{t ) := vj q (t)i dq (t) = V d (t)I d {t) + V q {t)I q (t). 


The instantaneous three-phase reactive power is defined as 


Q(t ) : = v dq(t ) 


1 

0 


idq(t ) = V d (t)I q (t) - V q (t)I d (t). 


Finally, the instantaneous three-phase (complex) apparent 
power is defined as 


S(t) :=P(t)+jQ(t). 


From the above definition, straightforward calculations 
together with standard trigonometric identities yield 


P(t) = 3 V(t)I(t) cos(0(f) — 99 (f)), 
Q(t ) = 3 V(t)I(t) sin(0(f) — 99 (f)). 


It follows that whenever v a b c and i a b c given in (2.4) pos¬ 
sess constant amplitudes, as well as the same frequency, 
i.e., 9 = ip, all quantities P, Q and S are constant. More¬ 
over, then the given definitions of power are in accordance 
with the conventional definitions of power in a symmetric 
steady-state [391141 [371. For further information on defini¬ 
tions and physical interpretations of instantaneous power, 
also under asymmetric conditions, the reader is referred to 

1133117311331 firm 131 H3T]. 


Since this work is mainly concerned with dynamics of 
generation units, all powers are expressed in “Generator 
Convention” [33 Chapter 2]. That is, delivered active 
power is positive, while absorbed active power is negative. 
Furthermore, capacitive reactive power is counted posi¬ 
tively and inductive reactive power is counted negatively. 


2.f. Algebraic graph theory 

An undirected graph of order n is a tuple Q := (V,£), 
where V := {m,..., n n } is the set of nodes and £ C [V] 2 , 
£ := {ei,..., e m }, is the set of undirected edges, i.e., the 
elements of £ are subsets of V that contain two elements. 
In the case of multi-agent systems, each node in the graph 
typically represents an individual agent. For the purpose 
of the present work, an agent represents a DG or storage 
unit, respectively a load. The Z-tli edge connecting nodes i 
and k is denoted by e; = {i, k}. By associating an arbitrary 
ordering to the edges, the node-edge incidence matrix B £ 
Rl v lx| £ l is defined element wise as bu = 1, if node i is 
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the source of the l -th edge ej, bn = — 1, if z is the sink 
of ei and bu = 0 otherwise. For further information on 
graph theory, the reader is referred to, e.g., [25] [35] and 
references therein. 

3. The microgrid concept 

Microgrids have attracted a wide interest in different 
research and application communities over the last decade 
[351 1351 1511 133] . However, the term “microgrid” is not 
uniformly defined in the literature [5T] [35] [33] [55] [T31 [37] 
133] . Based on [33] [35] [33] , the following definition of an 
AC microgrid is employed in this survey paper. 

Definition 3.1. A n AC electrical network is said to be an 
AC microgrid if it satisfies the following conditions. 

1. It is a connected subset of the LV or MV distribution 
system of an AC electrical power system. 

2. It possesses a single point of connection to the re¬ 
maining electrical power system. This point of con¬ 
nection is called point of common coupling (PCC). 

3. It gathers a combination of generation units, loads 
and energy storage elements. 

4. It possesses enough generation and storage capacity 
to supply most of its loads autonomously during at 
least some period of time. 

5. It can be operated either connected to the remain¬ 
ing electrical network or as an independent island 
network. The first operation mode is called grid- 
connected mode and the second operation mode is 
called islanded, stand-alone or autonomous mode. 

6. In grid-connected mode, it behaves as a single con¬ 
trollable generator or load from the viewpoint of the 
remaining electrical system. 

7. In islanded mode, frequency, voltage and power can 
be actively controlled within the microgrid. 

According to Definition |3.1[ the main components in a 
microgrid are DG units, loads and energy storage elements. 
Typical DG units in microgrids are renewable DG units, 
such as photovoltaic (PV) units, wind turbines, fuel cells 
(FCs), as well as microturbines or reciprocating engines 
in combination with SGs. The latter two can either be 
powered with biofuels or fossil fuels E3U3Z3- 

Typical loads in a microgrid are residential, commercial 
and industrial loads [M, 13 1531 . It is also foreseen to 
categorize the loads in a microgrid with respect to their 
priorities, e.g., critical and non-critical loads. This enables 
load shedding as a possible operation option in islanded 

mode [5X1153] . 

Finally, storage elements play a key-role in microgrid 
operation [53] [37]. They are especially useful in balancing 
the power fluctuations of intermittent renewable sources 
and, hence, to contribute to network control. Possible 
storage elements are, e.g., batteries, flywheels or superca¬ 
pacitors. The combination of renewable DGs and storage 


Main electrical network 



Figure 3: Schematic representation of a microgrid. The microgrid 
is composed of several DG units, loads and storage devices. The 
DG units are inverter-interfaced photovoltaic (PV), fuel cell (FC) 
and wind power plants. In addition, a power generation unit is con¬ 
nected to the network via a synchronous generator (SG). The point 
of connection of the microgrid to the main network is called point of 
common coupling (PCC). 


elements is also an important assumption for the inverter 
models derived in this paper. An illustration of an exem¬ 
plary microgrid is given in Fig. [3] 

Most of the named DG and storage units are either 
DC sources (PV, FC, batteries) or are often operated at 
variable or high-speed frequency (wind turbines, microtur¬ 
bines, flywheels). Therefore, they have to be connected to 
an AC network via AC or DC/AC inverters [40!, 101] . For 
ease of notation, such devices are simply called “invert¬ 
ers” in the following. Overviews on existing test-sites and 
experimental microgrids around the globe are provided in 
the survey papers [35] [5] [53] [53] [33] . 

Remark 3.2. While not comprised in Definition [XT] true 
island power systems are sometimes also called microgrids 
in the literature [35]. This can be justified by the fact 
that island power systems operating with a large share of 
renewable energy sources face similar technical challenges 
as microgrids. Nevertheless, an island power system differs 
from a microgrid in that it cannot be frequently connected 
to and disconnected from a larger electrical network m- 


Remark 3.3. Microgrids can also be implemented as DC 
systems [55][S3ED]- Definition 3.1 can easily be adapted to 
this scenario by removing the property “frequency control” 
in point [7] Recent reviews of the main differences and 
challenges for AC and DC microgrids are given in [ 1071 
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4. Modeling of inverter-based microgrids 

As outlined in the previous section, an AC microgrid 
is a spatially small power system the main components of 
which are renewable DG units, loads and energy storage el¬ 
ements interconnected through a network of AC transmis¬ 
sion lines and transformers. Furthermore, most renewable 
DG and storage units are interfaced to the network via 
inverters. As a consequence, fundamental network con¬ 
trol actions, such as frequency or voltage control, have 
to be performed by inverter-interfaced units. This fact 
represents a fundamental difference to the operation of 
conventional power systems, where mainly SG units are 
responsible for network control. Therefore and since the 
modeling of SGs is a well-covered topic in the literature 
[551 [5] HIJI, we focus in the following on microgrids with 
purely inverter-interfaced DG and storage units. Also, it 
is straightforward to incorporate SG-based units into the 
microgrid model presented hereafter. 

In line with these considerations, an inverter-based mi¬ 
crogrid can be represented by an undirected graph Q = 
(A/",£), in which—similarly to [Mj— nodes represent volt¬ 
age buses, edges represent dynamic power lines and the 
topology of the network is fully described by the incidence 
matrix B. The set of neighbors of node i is denoted by 
Mi and contains all k for which ej = {i,k} £ E. Please 
see Section |2.4| for a brief introduction to algebraic graph 
theory. In the present case, we further assume that the 
set of nodes M can be partitioned into two subsets Mi and 
Mr, associated to inverter and load nodes respectively. We 
next proceed as follows. First, we provide a description of 
the basic functionality and common operation modes of 
inverters that are instrumental for the modeling. There¬ 
after, we present models of inverters—depending on their 
mode of operation—loads, power lines and transformers. 
The section is concluded by combining the individual mod¬ 
els to an overall representation of an inverter-based micro¬ 
grid. To enhance readability, the subindex i, preceded by a 
comma when necessary, denotes in the sequel the elements 
corresponding to the *-th subsystem. 

4-1- Basic functionality and common operation modes of 
inverters 

Recall that inverters are key components of microgrids. 
Therefore, this section is dedicated to the model derivation 
of an inverter in a microgrid. The basic functionality of 
an inverter is illustrated in Fig. [4] The main elements of 
inverters are power semiconductor devices [SUES]. An ex¬ 
emplary basic hardware topology of the electric circuit of a 
two-level three-phase inverter constructed with insulated- 
gate bipolar transistors (IGBTs) and antiparallel diodes 
is shown in Fig. [5] The conversion process from DC to 
AC is usually achieved by adjusting the on- and off-times 
of the transistors. These on- and off-time sequences are 
typically determined via a modulation technique, such as 
pulse-width-modulation ED Eg. To improve the quality 
of the AC waveform, e.g., to reduce the harmonics, the 




Inverter 


Figure 4: Schematic representation of a DC-AC voltage conversion 
by a DC/AC inverter. The DC signal ‘ R>o —>• R on the left side 

is converted into an AC signal via power semiconductor devices. The 
generated AC signal uac,i : R>o R 3 a t the output of the power 
electronics is not sinusoidal. Therefore, an LC filter is connected in 
series with the power electronics to obtain an approximately sinu¬ 
soidal ouput voltage v\q 2 : R>o R 3 with low harmonic content. 



Figure 5: Typical circuit of a two-level three-phase inverter with 
LC output filter to convert a DC into a three-phase AC voltage. 
The inverter is constructed with insulated-gate bipolar transistors 
(IGBTs) and antiparallel diodes. The DC voltage is denoted by roc • 
R >0 —> R, the three-phase AC voltage generated by the inverter 
with v a b c : R> 0 —> R 3 , v a ^ c = col(u a , v^, v c ) and the three-phase 
grid-side AC voltage by vq : R>o R 3 > V G = col {vG a v Gbt v G c ) m 
The components of the output filter are an inductance Lf E R>o, a 
capacitance Cf £ R>o and two resistances Rf £ R>o, respectively 
Rf 2 £ R>o- Typically, the resistance R g £ R>o and the inductance 
L g £ M>o represent a transformer or an output impedance. At 
the open connectors denoted by “ 0 “ the circuit can be grounded if 
desired. 

generated AC signal is typically processed through a low- 
pass filter constructed with LC(L) elements. Further in¬ 
formation on the hardware design of inverters and related 
controls is given in [8T11751 112] . 

In microgrids, two main operation modes for invert¬ 
ers can be distinguished 11071 IBH : grid-forming and grid¬ 
feeding mode. The latter is sometimes also called grid¬ 
following mode [55j or PQ control [65], whereas the first 
is also referred to as voltage source inverter (VSI) con¬ 
trol [65]. The main characteristics of these two different 
operation modes are as follows [Ml j55J 1U7 , IMj . 

1. Grid-forming mode (also: VSI control). 

The inverter is controlled in such way that its out¬ 
put voltage can be specified by the designer. This 
is typically achieved via a cascaded control scheme 
consisting of an inner current control and an outer 
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voltage control as shown in Fig. 6a based on 


The feedback signal of the current control loop is 
the current through the filter inductance, while the 
feedback signal of the voltage control loop is the in¬ 
verter output voltage v a bc ■ R>o —i► R 3 - The inner 
loop of the control cascade is not necessary to con¬ 
trol the output voltage of the inverter and can hence 
also be omitted. Nevertheless, it is often included to 
improve the control performance. 

2. Grid-feeding mode (also: grid-following mode, PQ 
control). 

The inverter is operated as power source, i.e., it pro¬ 
vides a pre-specified amount of active and reactive 
power to the grid. The active and reactive power set- 
points are typically provided by a higher-level con¬ 
trol or energy management system, see mmm- 
Also in this case, a cascaded control scheme is usu¬ 
ally implemented to achieve the desired closed-loop 
behavior of the inverter, as illustrated in Fig. [6b] As 
in the case of a grid-forming inverter, the inner con¬ 
trol loop is a current control the feedback signal of 
which is the current through the filter inductance. 
However, the outer control loop is not a voltage, but 
rather a power (or, sometimes, a current) control. 
The feedback signals of the power control are the 
active and reactive power provided by the inverter. 


In both aforementioned operation modes, the current 
and voltage control loops are, in general, designed with 
the objectives of rejecting high frequency disturbances, en¬ 
hancing the damping of the output LC(L) filter and pro¬ 
viding harmonic compensation mZ> | T2, 73J HO]. Further¬ 
more, nowadays, most inverter-based DG units, such as 
PV or wind plants, are operated in grid-feeding mode |84j . 
However, grid-forming units are essential components in 
AC power systems, since they are responsible for frequency 
and voltage regulation in the network. Therefore, in mi¬ 
crogrids with a large share of renewable inverter-based DG 
units, grid-forming capabilities often also have to be pro¬ 
vided by inverter-interfaced sources [55., 13B]. 


Remark 4.1. Some authors [1071 [831 also introduce a third 
operation mode for inverters called grid-supporting mode. 
Nevertheless, this last category is not necessary to classify 
typical operation modes of inverters in microgrids in the 
context of this work, since grid-supporting inverters are 
grid-forming inverters equipped with an additional outer 
control-loop to determine the reference output voltage. 
Such outer control-loops are, e.g., the usual droop controls 
H7H3H- Therefore, the term “grid-supporting inverter” is 
not used in the following. 


Remark 4.2. In addition to the two control schemes in¬ 
troduced above, there also exist other approaches to op¬ 
erate inverters in microgrid applications. For example, 
[5J [106; 111] propose to design the inverter control based 
on the model of an SG with the aim of making the in¬ 
verter mimic as closely as possible the behavior of an SG. 


However, to the best of the authors’ knowledge, these ap¬ 
proaches are not as commonly used as the control schemes 
shown in Fig. |6a| and Fig. 

4-2. Modeling of grid-forming inverters 

A suitable model of a grid-forming inverter for the pur¬ 
pose of control design and stability analysis of microgrids 
is derived. There are many control schemes available to 
operate an inverter in grid-forming mode, such as PI con¬ 
trol in (^-coordinates [81] , proportional resonant control 
[551 TOO] or repetitive control 1081 [50] among others. An 
overview of the most common control schemes with an 
emphasis on H a 0 repetitive control is given in [112] . For 
a comparison of different control schemes, the reader is 
referred to [53]. The assumption below is key for the sub¬ 
sequent model derivation. 

Assumption 4.3. Whenever an inverter operated in grid¬ 
forming mode connects a fluctuating renewable generation 
source, it is equipped with a fast-reacting storage. 

Assumption |4.3| implies that the inverter can increase 
and decrease its power output within a certain range. This 
is necessary if the inverter should be capable of providing 
a fully controllable voltage also when interfacing a fluc¬ 
tuating renewable DG unit to the network. Furthermore, 
since the storage element is assumed to be fast-reacting, 
the DC-side dynamics can be neglected in the model. The 
capacity of the required DC storage element depends on 
the specific source at hand. Generally, the standard ca¬ 
pacitive elements of an inverter don’t provide sufficient 
energy storage capacity and an additional storage compo¬ 
nent, e.g., a battery or flywheel, is required if an inverter 
is onerated in grid-forming mode 1101 See m for a sur¬ 
vey of energy storage technologies in the context of power 
electronic systems and renewable energy sources. 

Due to the large variety of available control schemes, 
it is difficult to determine a standard closed-loop model of 
an inverter operated in grid-forming mode together with 
its inner control and output filter. Therefore, the ap¬ 
proach taken in this work is to represent such a system 
as a generic dynamical system. Note that the operation 
of the IGBTs of an inverter occurs typically at very high 
switching frequencies (2-20 kHz) compared to the network 
frequency (45-65 Hz). It is therefore common practice 
[53 0E EH CZ2CES] to model an inverter in network studies 
with continuous dynamics by using the so-called averaged 
switch modeling technique EH HE], be., by averaging the 
internal inverter voltage and current over a suitably chosen 
time interval such as one switching period. 

It is convenient to partition the set J\fj into two subsets, 
i.e., A/j = A/i U AT, such that AT contains all nodes asso¬ 
ciated to grid-forming inverters and AT contains those as¬ 
sociated to grid-feeding inverters. Consider a grid-forming 
inverter located at the i-th node of a given microgrid, i.e., 
i e A/i- Denote its three-phase symmetric output voltage 


(>b 
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(a) Schematic representation of an inverter operated in grid-forming 
mode based on ED- Bold lines represent electrical connections, while 
dashed lines represent signal connections. The current through the 
filter inductance is denoted by if, a bc : R>o R 3 and the inverter 
output voltage by v a b c : M>o —>• R 3 . Both quantities are fed back 
to a cascaded control consisting of an outer voltage and an inner 
current control. The reference signal u re f : R>q —> R 3 for the voltage 
controller is set by the designer, respectively a higher-level control. 
The IGBTs of the inverter are then controlled via signals generated 
by a modulator. The control structure can also be reduced to a pure 
voltage control. 



(b) Schematic representation of an inverter operated in grid-feeding 
mode based on ED Bold lines represent electrical connections, while 
dashed lines represent signal connections. As in Fig. |6a| the current 
through the filter inductance is denoted by if a b c : R>q — > R 3 and the 
inverter output voltage by u a j, c : R>q —> R^. In grid-feeding mode, 
both quantities are fed back to a cascaded control consisting of an 
outer power and an inner current controller. The reference active 
and reactive powers P re f E R, respectively Q re f E R, are set by the 
designer or a higher-level control. 


Figure 6: Schematic representation of the two main operation modes of inverters in microgrid applications: grid-forming and grid-feeding 
mode. 


by v a bc,i ■ R>o — X R 3 with phase angle a* : R>o —X S and 
amplitude R> 0 —x R> 0 , Le., 


Vabc.i — \l ^ ^ 


sin(a,;) 
sin(cq — 1 7 r) 
sin(aj + §7r) 


Furthermore, denote by w* := on the frequency of the volt¬ 
age v a bc,i • Denote the state signal of the inverter with its 
inner control and output filter by 27 : R>o —X R m , its in¬ 
put signal by u re f y : R>o — X R 3 and its interconnection 


, see Fig. 


6a 


port signals by v a bc,i and i a bc,i ■ R>o - 
Let ft : R m x I 3 x I 3 -> K m and hi : R m x I 3 -) R 3 de¬ 
note continuously differentiable functions and i\ denote a 
nonnegative real constant. Then, the closed-loop inverter 
dynamics with inner control and output filter can be rep¬ 
resented in a generic manner as 


^ref,iDa&c,*)) 
'Oabc,i — h,; (Xi , V re f_ j ), 


(4.1) 


where the positive real constant 7 denotes the time-drift 
due to the clock drift of the processor used to operate 
the inverter, see [9TJ for further details. Note that i a bc,i 
represents a disturbance for the inner control system. 

One key objective in microgrid applications is to design 
suitable higher-level controls to provide a reference voltage 
u re f,i for the system (4.1) m- Within the hierarchical 


control scheme discussed, e.g ., in mm this next higher- 
control level corresponds to the primary control layer of a 
microgrid. Let Zi : R>o —X R p denote the state signal of 
this higher-level control system, Ui : R>o —X R 9 its input 
signal and v re fi its output signal. Furthermore, let gi : 
]R P x ]R 9 —» R p and Wi : R p x R 9 —x R 3 be continuously 
differentiable functions. Then, the outer control system of 
the inverter can be described by 


7 ik = 9i(zi,Ui), 
v re { ,i = Wi{Zi,Ui). 


(4.2) 


Combining (4.1) and (4.2) yields the overall inverter dy¬ 
namics for the i-th node, i £ ATi , 


7 iii = gi{zi,Ui), 

Vabc,i = hi(Xi,Wi(Zi,Ui)). 


(4.3) 


4-3. Modeling of grid-feeding inverters and loads 

As discussed in Section |4.1| grid-feeding inverters are 
typically operated as current or power sources. In order 
to achieve such behavior, the control methods employed to 
design the inner control loops of grid-forming inverters (see 
Section 4.2) can equivalently be applied to operate invert¬ 
ers in grid-feeding mode. The current or power reference 
values are typically provided by a higher-level control, e.g., 
a maximum power point tracker (MPPT) |84| . 
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We define the set Aft := A /2 U AFr that contains the 
nodes associated to grid-feeding inverters and loads. As 


done for the model of a grid-forming inverter in (4.3), let 
xej : R>o —X R r denote the state signal, v a bc,i ■ R>o —► R 3 
and i a bc,i '■ R>o —X R 3 denote the interconnection port 


signals, : R r x R d -) R r and h^i : R r x R J -> R d 
denote continuously differentiable functions and Ki denote 
a nonnegative real constant. We assume then a generic 
dynamic model of the form 


KiX£,i — fc,i {x£,i > iabc,i) ; 
Vabc,i — iabc,i )? 


(4.4) 


for any node i £ Aft- In addition to grid-feeding invert¬ 


ers, the model (4.4) can equivalently represent impedance 


( e.g ., R parallel to L), current- or power-controlled loads. 
Furthermore, a large variety of other load behaviors can 


be modeled by (4.4). We refer the reader to [59, 103j for 
further details on load modeling. 

4-4- Modeling of power lines and transformers 

The main purpose of the present paper is to provide a 
structured modeling procedure for microgrids. For ease of 
presentation, we make the following assumption. 

Assumption 4.4. All power lines and transformers can 
be represented by symmetric three-phase RL elements. 

In light of Assumption |4.4| and to ease presentation, 
we solely use the term power lines to refer to the network 
interconnections in the following sections. Also, note that 
it is straightforward to extent the modeling approach pre¬ 
sented hereafter to more detailed power line or transformer 
models, as well as to DG units interfaced to the network 
via SGs. 

Recall that the topology of a microgrid can be con¬ 
veniently described by an undirected graph Q = (A/",£), 
where £ denotes the set of power lines interconnecting 
the different network nodes i ~ A f. We associate an ar¬ 
bitrary ordering to the power lines e; ~ £. Likewise, we 
assign to each power line ei € £ a three-phase line current 
With Assumption 


lL,l 


^>o 


4.4 


the power line 
£ [A f] 2 is sym- 


ei £ £ connecting a pair of nodes {*, k} 
metric, i.e., each phase of the power line ej is composed 
of a constant ohmic resistance Ri G R>o in series with a 
constant inductance Li G R>o and Ry as well as Li have 
the same value for each phase. Furthermore, the voltage 
drop across the line is given by 

Vl,1 * = Vabc.i 'Vabc.k- 

We denote the state of the Z-th line by xl,i ■= ih,i and its 
interconnection port variables by vl,i and Then, the 
model of the Z-tli line ei G £ is given by 


Li±l,i = -Rix l ,i + vl,i, 
il,i = x L ,i- 


(4.5) 


For a compact derivation of the network dynamics, it is 
convenient to define the aggregated nodal voltages and 
currents 

Vabc ■= col [v a bc,i) £ R 31 ^ 1 , i a b c ■= col (i a bc,i) £ R 31 ^ 1 , 
the aggregated line voltages and currents 

v L := col (v L ,i) £ R 3|£| , %l := col(i L ,i) £ R 3|£| , 
as well as the matrices 

L := diag(Lz) G R |£:|x|£:| , R := diag(Rz) £ R |£r|x|£:| . 

Then, the three-phase interconnection laws can be ob¬ 
tained by following the approach used in |34l , where Kirch- 
hoff’s current and voltage laws (KCL and KVL) are ex¬ 
pressed in relation to the node-edge incidence matrix B, 
i.e., 


^abc — B I 3 Z/,, 


B 


> hv a bc = VL- 


(4.6) 


Hence, by combining (4.5) with (4.6) the dynamical system 


representing the network is given by 

L <g> 1 3 ±L = -R <8 I 3 XL + B 1 
^abc — B 8 T-3XL- 


J ^-3Vabc: 


(4.7) 


We next transform the model (4.7) into dq-coordinates 


by means of the transformation Td q introduced in (2.3). 
This coordinate transformation is instrumental for the model 
reduction carried out in Sectional Let 


:= mod27r (w com Z) £ S, 


(4.8) 


where the operatoi)^] mod 27I . is added to respect the topol¬ 
ogy of the torus. Applying the transformation Td q with 
transformation angle <j> to the signals v a i , c y and i a bc,i, i ~ 
A f, gives 


Vdq.i ’■— 'A c i q {4 i )v a bc,i — 


Vd,i 


v„ 


q 


j Idq.i • — r £dq((f')iabc,i — 


Id,i 


L qo. 


where the superscript ” : ” is introduced to denote signals 
in dq-coordinates with respect to the angle </>. This nota¬ 
tion is used in the subsequent section, where a reduced- 
order model of a microgrid is derived by using several dq- 
transformation angles. Furthermore, following standard 
notation in power systems, the constant </> = w com is re¬ 
ferred to as the rotational speed of the common reference 


frame. Likewise, the signal Xl,i in (4.51 becomes 


X L.dq.l •— r £dq(.4 > ')XL 1 abc,l — 


x Lt0 


Note that 


XL.dq.l — 'i'd q {.f£)XL,l T Tdqi.fyxL,l 


= 00 


-X 


L,q,l 


X 


L,d,l . 


Td q (4>)XL,l- 


3 The operator mod 27 r : R —X [0,27r), is defined as follows: y = 
rnorpTT {■>'} yields y = x—k2tr for some integer k, such that y S [0, 2tt) . 
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Hence, (|4.5|) reads in dg-coordinates as 


1 , .com 


r 

-1 

• "S 

_1 


LlXL,dq,l = L >1 

= ~RlXL,dq,l + LlUJ C 
iL,dq,l — %L,dq,l’ 


—X L ,q,l 


X 


L,d,l 


“t“ Xdq,l , 


(4.9) 


By defining the aggregated nodal voltages and currents in 
dg-coordinates 

Vdq := CO \{v d q,i) G R 2|A ^, idq ■= COl {i d q,i) G , 


Network 


L (g> l 2 XL,dq = (-R (8) I 2 + X)xL,dq +B T ® 1 2 V d q, 
^dq = I|dV| & Tdq((f))V a b c , 

iabc — 13 0 Tdq{$) X^^q 


Vabc,i 


^ abc,i 1abc,k 


Vabc,k 

i cj at 

7Vi = gi{zi,Ui), 

Viliii = fi(xi,Wi(zi,Ui),i a ^i), 

Vabc,i = hi{Xi,Wi(Zi,Ui)) 

. 

k ~ J\fi 

ftkXl^k — f/’,k\X£,ki'iabc,k)i 

Vabc,k = h£,k(x£fi, iabc,k) 

. 


Grid-forming inverters Grid-feeding inverters and loads 


the aggregated line voltage and currents in (^-coordinates Figure 7: Port-b as ed rep resentation of an inverter based-microgrid 

modeled by ( |4.11| ), ( |4.12| t. 

VL,dq ■= col (v L>dq ,i) G K 2|£| , x L>dq := col {i L ,dq,i) G K 2|f| , 


as well as the matrix 


X := diag 



0 

1 



g r2|£|x2|£|^ 


(4.7) becomes in dq- coordinates 


L 8 l 2 x L ,dq = (-R ® h + X)x L , dq + B t <8 1 2 Vdq, 

idq — H (8 I 2 %L,dq- 

(4.10) 


4-5. Overall model 

By defining the state vectors z G RP^I, x G ]R TO I- A/ " 1 I ) 
Xu G the input u G the matrices 

T : = diag( 7 *) G R |M|x|M| , v ;= diag (*,.) £ RlMMAil, 
dT: = diag(K i )GK |A/ ' 2|x|A/ ' 21 , 


and combining (4.31, (4.4) and (4.10), the overall microgrid 
model (see Fig. |7j) is given by the differential equations 


[r® I p |jvi|] z = g(z, u), 

[vr 8 I m |Ad|] * = f(x,w{z,u),i abc ), 

\X (8 Ir |A/11 ] f£ J 'labc') > 

[L 8 I 2 ] x L>dq = (-R 8 I 2 + X)x LAq + B t 8I 2 v dq , 


(4.11) 


together with the algebraic relations 


^ abc — H 8 I 3 Xdj — 13 8 T dq ((j)) XL^ dq , 

Vabc,i — hi{Xi,Wi(Zi,Ui)), 

Vabc,k (X&,k 7 labc,k) 5 

Vdq = Vl 8 T dq ((f>)v abc , * ~ M, k ~ TV). 


(4.12) 


5. To phasors and voltage sources via time-scale 
separation 

For the purpose of deriving an interconnected network 
model suitable for network control design and stability 
analysis, it is customary to make the following assump¬ 


tions on (4.11), (4.12), where e stands for a generic small 


positive real constant. 


Assumption 5.1. Vi < e in (4.3), i ~ TV). Therefore, 


ii{t ) = 0 m for all t > 0. Furthermore, v abc d = Wi(zi,Ui), 
i ~ TV). 


Assumption 5.2. K k < e in (4.4), k ~ TV). Therefore, 


xe,k(t) = 0 r for all t > 0. Furthermore, the power balance 
at each node k G TV) can be described by a ZIP model m 


i.e., 


Pk{v d q,k,idq,k ) = ^ {ap,k\\v dq ,k\\ 2 + bp,k\\vdq,k\\2 + Cp,fc) 

Qk(x d q,k, l d q,k) — || V d q,k || 2 “b || Vdq,k || 2 V CQ ik ) 

’■= Q*k(\\Vdq,kh), 


where ap jfc , b P ^ k , cp jfc , aQ, fc , b Qik and CQ ik are real con¬ 
stants and P k {v dq ,k,i dq ,k) and Q k (vd q ,k,idq,k ) are calcu¬ 
lated as given in Definition\2.11\ 


lumption 5.3. L < el|£| in (4.10). 


Therefore, 


Assumption |5.1| is equivalent to the assumption that 
the inner current and voltage controllers track the voltage 
and current references instantaneously and exactly. Usu¬ 
ally, the current and voltage controllers in (|4.1|) (see also 


Fig. 6a) are designed such that the resulting closed-loop 
system (4.1) has a very large bandwidth compared to the 
control system located at the next higher control level rep¬ 
resented by (4.2) |65j 2Ql 74]. If this time-scale separation 
is followed in the design of the system (4.3), the first part 
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of Assumption |5. l| can be mathematically formalized by in¬ 
voking singular perturbation theory [5BJ Chapter 11], [55] . 
The second part of Assumption |5.1| expresses the fact that 


the inner control system (4.1) is assumed to track the ref¬ 


erence v Te f'i = u>i(zi,Ui ) exactly, independently of the dis¬ 
turbance i a bc,i- Typical values for the bandwidth of (4.11 
reported in [7i m] are in the range of 400 — 600 Hz, while 
those of ( |4.2[ ) are in the range of 2 — 10 Hz. 

Assumption |5 . 2| implies that the dynamics of loads and 
grid-feeding units can be neglected. This assumption is 
also frequently employed in microgrid and power system 
stability studies, where loads are often modeled as either 
constant impedance (Z), constant current (I) or constant 
power loads (P) or a combination of them (ZIP) [5511103] , 
Similarly, grid-feeding units with positive active power in¬ 
jection are represented by setting ap } k = ciQ k = 0 and 
bp,k or c P,k to negative values. The values for bqj. and 
CQ,k should be chosen in dependency of the reactive power 
contribution of the unit. 


Assumption 5.3 is standard in power system analysis 
|52| [321030152103- The usual justification of Assump¬ 
tion |5.3| is that the line dynamics evolve on a much faster 
time-scale than the dynamics of the generation sources. In 
the present case, Assumption |5.3| is justified whenever As¬ 
sumption 5.1 is employed, since the line dynamics (4.101 


are typically at least as fast as those of the internal inverter 
controls (4.1), see, e.g., [81]- Again, Assumption 5.3 


can 


be mathematically formalized by invoking singular pertur¬ 
bation arguments [55] Chapter 11], [55] . 

Under Assumption |5. 1 [ the model of each grid-forming 


inverter (4.3) reduces to 


'7i~i — 9i( z ii u i) i 
^abc,i — U’i)i ^ ^ M 1 • 


(5.1) 


The model (5.1) represents the inverter as an AC voltage 


source, the amplitude and frequency of which can be de¬ 


fined by the designer. The system (5.1) is a very commonly 


used model of a grid-forming inverter in microgrid control 
design and analysis [55] @0] [55] [50]. 


Furthermore, often a particular structure of (5.1) is 


used in the literature [Ml ESI M HH El [76]. As discussed 
in Section |2.1[ a symmetric three-phase voltage can be 
completely described by its phase angle and its amplitude. 
In addition, it is usually preferred to control the frequency 
of the inverter output voltage, instead of the phase angle. 
Hence, a suitable model of the inverter at the i-th node is 
given by I5D] !HS] 


"iidu =Uli = U f, 

Vi = uj, 

Vabc,i Vabc,i 7 U) i 


(5.2) 


where uf : R>o —> R and u( : R>o —> R are control signals. 

Usually, it is also assumed that the active and reactive 
power output of the inverter is measured and processed 


^ref 


^ref 



Figure 8: Simplified representation of an inverter operated in grid¬ 
forming mode as ideal controllable voltage source. Bold lines repre¬ 
sent electrical connections, while dashed lines represent signal con¬ 
nections. 


through a filter to obtain the power components corre¬ 
sponding to the fundamental frequency [SHIM] EH 


(5.3) 


Here, Pj and Qi are the active and reactive power injec¬ 
tions of the inverter, P™ : R> 0 —► R and Q[" : R>o —> R 
their measured values and rp i £ R>o is the time constant 
of the low pass filter. 



Note that whenever the particular form (5.2), (5.3) of 


(4.2) is considered and the measured and filtered power 


signals are used as feedback signals in the controls u". re¬ 
spectively uY, then the bandwidth of the overall control 
system is limited by the bandwidth of the measurement 


filter. Hence, if Tp i j/,, then Assumption 5.1 is justified. 


With Assumption 5.2 
algebraic relation 


(4.4) can be represented by the 


Pk{Vdq,kPdq,k) — Pk (|| Vdq,k H2), 
Qk{ydq,ki'id.q,k) — Qk (I! Pdq,k || 2)7 


k ~ Mi. 


(5.4) 


Finally, under Assumption |5.3| the network model (4.10) 
is also static and given by 


i dq = B® I 2 (i? <8> I 2 — A”) -1 B r < 


1 t-2Vdq- 


(5.5) 


The reduced-order model (5.2) - (5.51 is still rather 


complex to handle, as the variables of (5.2) are expressed 


in a6c-coordinates, while those of (5.4) and (5.5) are ex¬ 


pressed in common dg-coordinates. Therefore, a more 


compact representation of (5.2) - (5.5) is derived in the 
following. To this end, it is convenient to recall that ctj 
is the angle of the voltage at the i-th node with initial 
condition a 0i , i ~ M and to define 


Si := op; + [ {o-i — L0 cora )dT e S, i ~ M. 
Jo 


(5.6) 
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Vabc,i 



Figure 9: Illustration of the different coordinate frames used to 
derive the model of an electrical network given in The signal 

v abc i '■ M> 0 —> M 3 denotes the three-phase voltage at the i-th bus 
with phase angle a* : M>o —> § and amplitude Vi : R>o —> M>o, 
i J\f. The mappings Tdq and T$ are given in J2T3J, respectively 
The angle 5i : M>o —> S is defined in ( |5.6| T ' Note that, by 
construction, — Si = mod 27 r (ui com t) = <f>, where the real constant 
cjcom denotes the speed of the common dg-reference frame. 


Let w : M> 0 —> S and consider the mapping T$ : § —> 


T s (w) := 


cos(n7) sin(n7) 
— sin(tu) cos(t*7) 


(5.7) 


Note that, with <5, defined in ( |5.6| ), 

OLi — Si = mod 2 ,r (w com t) = <j>, i ~ TV, 
and that straightforward algebraic manipulations yield 
Tdq(<j>) = T 5 (6i)T dq (ai). 


Note that the reactances Xi = Liui com are calculated at 
the frequency w com , which, under the made assumptions, 
should be chosen as the (constant) synchronous frequency 
of the network—denoted by uj s £ M in the following] Typ¬ 
ically, lo s £ 27r[45, 65] rad/s. 


Remark 5.4. The form (5.12) is a very popular represen¬ 
tation and these complex quantities are often denoted as 
phasors (5l Hl2|. Furthermore, by using Euler’s formula 
HZ], p~12| ) can also be rewritten in polar form. Note, 
however, that, unlike, e.g., ©uni, other authors define a 
phasor as a complex sinusoidal quantity with a constant 
frequency m- 


Define the admittance matrix of the electrical network 


by 

y~B{R + jX)~ 1 B r eC |Af|x|Af| 

(5.14) 

and 

:=K(y«), B u :=Z{y u ), 

Vik •— G T jB^k .— Viki i 7 ^ k. 


Moreover, it follows immediately that 



1° 

y* = \ 


if nodes i and k are 
not connected 
if nodes i and k are 
connected by line l 


Hence, by construction, 

^dq.i — r Rdq{^')'^abc,i — 'Rs( < &i)'Vdq(,OLi) r Vabc,i — Rs^&i'jVi ^ ; 

(5- 8 ) 

which makes it convenient to define 


and 

Gn + jBa = {Ri + jXi) 1 , 

l~£i 


where denotes the set of edges associated to node i. 
Inserting (5.10) and ( |5.11[ ) into ( |5.13[ ) yields 


Vdq,i ■= 


Vd,i 

— 1/ 

'o' 

v qA 

— Vi 

1 




(5.9) 


The variables v dqd are referred to as local dg-coordinates 
of v a bc,i in the following. The relation between v a bc,i, v dqd 
and v dqt i is illustrated in Fig. [9] 

It is convenient to represent ( |5.8[ ) in the complex plane 

Vqd,i := Vq,i T jV d ,i = (cos(d^) T j sin(dj)) V qd i = e^ l V qdd) 

(5.10) 

where V qdt i = V q ^ + jV d ,i, i ~ Af. Equivalently, let 


Iqd,i •— Iq.i T j^d,i — 1 I qd i (5.11) 


Iq d = diag (e jSi ) JMiag (. e jSi ) V qd . 


(5.15) 


Recall that V qd and I qd defined in ( |5.12[ ) are expressed in 
local dg-coordinates. By making use of (5.9) and (5.14), 
(5.151 can be written component-wise as 

Iqd,i — Iqd "F 3 Rl.i , 

Iqd = GaVi 'y ] {Gik cos(difc) T B^ sin(djfc)) V&, 

kr^Mi 

Id,i — B^Vi y ' (Bik COs(Sik) sin(d^fc)) Vfc, 

(5.16) 


and define 

Vqd ~Vq + jV d £ CM, iqd := Iq + jl d £ C'^l , 
Vqd ~V q + jV d £ C M , I qd ~ I q + jl d £ C 1 ^' 1 . 


(5.12) 


Then, with X := diag(Xj) = diag(L ; w com ) £ IRl £ l x l £ l, 
we can rewrite (5.5) as 


I qd = B(R + jX)~ 1 B T V qd . 


(5.13) 


i ~ Al, where, for ease of notation, angle differences are 
written as Sik := Si — 5k- Furthermore, from Definitional 1| 


4 Under the made assumptions, ([ 5 T 5 J is the equilibrium of the 
’’fast” line dynamics ( |4.9| ) |561 Chapter 11]. Hence, in order for the 
currents i^ q and voltages v^q to be constant in steady-state, a; com 
has to be chosen identically to the synchronous steady-state network 
frequency. 
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together with (5.9) and (5.16), the power flows in the net- 


Network 


work are given by 

Pi = Vil q< i = 

GuVf - 22 ( G ik cos (Sik) + B ik sm{8ik)) V k V i: 

k~Ni 

Qi — 

= —BaVi + 22 (Bik cos (Sik) - Gik sin (5 ik )) V k Vi. 

k~Jfi 

(5.17) 


* h! 


P t =GuV? - J2 (GA-cosfe) + B ik sin(<5is,)) V k V it 

k~Afi 

Qi = k»ii V-~ + ^ ( ( Bik. COs(<5 ik) Gi k V k Vi 

k~Afi 


{Si, Vi) {Pi, Qi) 


( P k ,Qk ) 


{S k ,V k ) 


The equations (5.17) are the standard power flow equa¬ 


tions used in most recent work on microgrid control design 
and stability analysis, e.g., [5511501 121175], 

Remark 5.5. Note that for any other choices of the trans¬ 
formation angle in local dg-coordinates V,i, p 0. This is 
usually the case when modeling SGs, since the angle of 
the internal machine electromagnetic force (EMF) is in 
general not known. Then, the equations (5.17) become 
slightly more involved, see [5] Chapter 9]. 


i ~ N\ 

7 A — uj- w com , 

li T Pi PY l = -P;r + Pi, 

Vi = uj, 

i,rp,QT = -QT + Qi 


k ~ Hi 


0 = P k - Tt (14), 
0 = Qk - Qt(V k ) 


Grid-forming inverters 


Grid-feeding inverters and loads 


Figure 10: Reduced microgrid model in standard representation 


with grid-forming inverters modeled by ( 5.18| ), as well as loads and 
grid-feeding inverters given by ( |5.4[ ) with | Vdq,k II 2 — Vk - The network 
is represented by the power now equations (|5.17[ ). This compact 
model is obtained from the model - i fruif by expressing the 

current and power flows in local dg-coordinates Vd q ,i and id q ,i, i ~ A/*, 
see ( |5.15| ) and Fig.[9j 


Furthermore, in local dg-coordinates, the particular in- Then the system (5.2) - (5.51 can be written equivalently 


verter model (5.2), (5.3), is given by 


by means of (5.4), (5.17), (5.18) as 


= Ui - a/ om = uf - w com , 

liT P ,P™ = -PT + Pi, 

V t = uY, 

7 iT Pi Q? = -Q? + Qi, 


(5.18) 


with V q d t i = Vi (see (5.9)) and Pi and Qi given by (5.17). 
Finally, recall (5.4) and note that 


H'Cdq.felb = ||V^d,fc||2 = HVgd.felb = V k , k £ A ft- 


This completes the reformulation of the model (5.2) - (5.5). 


The final overall microgrid model is given by (5.4), (5.17) 


T'\ r „ d , .com-ti 

Tdi = u — uj 1 |a/i | > 

TT pm = _pm + pi 

Vl ~ U ’ (5.19) 

TTQ m = -Q m + Q/, 

Vd = p e~ p h 

o m =Qt-Ql 

where the last 2|7V) algebraic equations correspond to the 
power balances at nodes k ~ Aft- 

This section has illustrated the main modeling steps 
and assumptions, which lead from the detailed microgrid 


(5.18) and shown in Fig. 10 This is the standard model model (4.11), (4.12) to the model (5.19), (5.17), respec- 


employed throughout the literature. 

The section is concluded by deriving a vector-based 


tively (5.4|, (5.17), (5.18). The model (5.19), (5.17) is 


formulation of the microgrid model (5.4), (5.17), (5.18). 
To this end, we define the vectors 


<5/ : 

=col(5i) G S |M| , 

Vf 

= col(Vi) 

u s : 

=col(uf) e R |MI , 

u v 

= col(wf) gK |M| , 

Pi : 

=col(Pi) G R |M| , 

Qi 

= col(Qi) G M |M| , 

Pt : 

=col(P fc ) G R |A/ 4 

Qi 

= col(Q fc ) 6l WI , 

Pt ■ 

=col(P fe *(V fc )) gm'^i, 

QI := co\(QUV k )) 


frequently used in the analysis and control design of mi¬ 
crogrids [M] (HU HU [Ml [Ml HI EH EU [Ml • Some of the men¬ 
tioned work is conducted under additional assumptions 
such as instantaneous power measurements j95j|2j[76], con¬ 
stant voltage amplitudes [951 [14))') S [89] or small phase angle 
differences [Ml 1521 155] , In addition, ideal clocks are usually 
assumed, i.e., T = Ii y, i. Furthermore, whenever constant 
impedance or constant current loads are assumed, the al¬ 


gebraic equations in (5.19), (5.17) can be eliminated by an 


with Pi, P k , Qi, and Q k given by (5.17), as well as the 
matrix 

T := diag( t P i ) £ R|Vi|x|A/\| 


appropriate network reduction. This process is commonly 
known as Kron reduction and frequently employed in mi¬ 
crogrid and power system studies. For further details on 
Kron reduction, the reader is referred to [Ml EH] . 
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6. Conclusions and topics of future research 

6.1. Summary 


The present survey paper has introduced the reader to 
the microgrid concept with the main focus of providing 
a detailed procedure for the model derivation of a three- 
phase inverter-based microgrid. In particular, it has been 
shown how—and under which assumptions—the microgrid 
models usually used in the literature can be obtained from 
a significantly more complex model derived from funda¬ 
mental physical laws. The assumptions invoked in the 
reduction process are often satisfied in standard applica¬ 
tions. Therefore, the reduced model represents a valid ap¬ 
proximation and may, hence, be useful for control design 
and system analysis. In addition, the employed model re¬ 
duction techniques can equivalently be applied to standard 
bulk power system models. 

Nevertheless, it is important to note that the model 
derived in the present paper neglects effects such as asym¬ 
metric operation, DC-side dynamics of DG units or line 
capacitances. These facts have to be kept in mind, when 
performing microgrid analysis based on the derived model 
and assessing the results. 

Also, it is worth mentioning that numerical simulation 
of the introduced microgrid models (4.111, (4.12), respec¬ 
tively (5.19), (5.17), requires careful selection of the nu¬ 
merical integration method to be employed. The main 
reason for this is that the model (4.11), (4.12) contains 
dynamics evolving at a wide range of time-scales, i.e., it is 
a stiff model [49j Chapter 7], [33] Chapter 8]. As is well- 
known, certain (standard) numerical integration methods 
will lead to numerical instability, when applied to stiff 
models—unless an extremely small step size is employed 
[321 Chapter 7], [331 Chapter 8]. On the contrary, in the 
reduced model (5.19), (5.17) the fast dynamics have been 
eliminated and replaced by their corresponding steady- 
state equations. Hence, this model is not stiff and simpler 
integration methods can be used. Very similar situations 
are usually encountered in simulation of large conventional 
power systems [521 Chapter 13], [311 . 


6.2. Future research 

To conclude this survey, we very briefly highlight some 
topics of future research. Following up on the discussion 
in the introduction, there are numerous challenges related 
to system and control theory in microgrid applications. 
To further motivate these, we briefly review control goals 
in microgrids. At the present, the following are consid¬ 
ered to be among the most relevant control objectives in 
microgrids [32[ [5TJ [301 SSI [551 EZ1 EI] : frequency stabil¬ 
ity, voltage stability, operational compatibility of inverter- 
interfaced and SG-interfaced units, desired power sharing 
in steady-state, seamless switching from grid-connected to 
islanded-mode and vice-versa, robustness with respect to 
uncertainties and optimal dispatch. 

Some of these problems have been addressed in recent 
work within the control community, e.g., frequency and 


voltage stability [251 SSI SOI 11021 1351 [31 [751 SH S31 132] , 
secondary control HOI S3 US CM S3] or optimal dispatch 
[321 EH [331- Compared to the model derived in the present 

paper, most of the aforementioned work is conducted un¬ 
der certain additional simplifying assumptions, such as 
constant voltage amplitudes, identical DG unit dynamics, 
as well as lossless or identical lines. 

On the modeling side, one direct extension of the pre¬ 
sented modeling framework is to investigate the suitabil¬ 
ity of the use of more refined modeling techniques such 
as dynamic phasors [381 HOI EH or symmetric components 
HUGH, to describe the dynamics of a microgrid in asym¬ 
metric operating conditions. In this context, one main 
challenge is to accurately consider these phenomena while 
maintaining an analytically tractable model. Another sub¬ 
ject of future research is the derivation of more detailed 
load models for microgrid applications. As in conventional 
power systems, accurate load modeling is a very impor¬ 
tant, but also very difficult task [5^11103] . The main rea¬ 
son for this is that there are typically many different kinds 
of loads connected within one power system or microgrid, 
see, e.g., [59l Chapter 7]. As a consequence, it is difficult 
to obtain suitable generically valid abstractions. 

In conclusion, there are many challenging open research 
questions regarding a reliable, safe and efficient opera¬ 
tion of microgrids. Therefore, the authors hope that the 
present survey on modeling of microgrids may serve as a 
stimulating base for a large variety of future research on 
both the theoretical and the application side. 
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